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Abstract 

We study the dynamics of bubble wall in five-dimensional Gauss-Bonnet gravity in the thin- 
wall approximation. The motion of the domain wall is determined by the generalized Israel 
junction condition. We find a solution where the wall of true vacuum bubble is oscillating, the 
breathing bubbles. We briefly comment on how this breathing vacuum bubble can affect the 
analysis of the string theory landscape. 
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1 Introduction 



In the theory with multiple meta-stable vacua, string theory for example, our universe is very 
likely to tunnel from other vacua with higher cosmological constants. The process is similar to 
the first order phase transition, where the bubble with different vacuum energy will nucleate from 
the original vacuum. When couple to the gravity, the vacuum with positive energy will expand 
exponentially. This leads to an interesting cosmology of eternal inflation [l]. 

In Einstein gravity, the geometry and dynamics of the bubble nucleation have been studied 
extensively, see [3] and references therein. The motion of the bubble wall, in the thin wall approxi- 
mation is determined by the Israel junction condition [2j. There are various kinds of solutions, but 
they fall into two main categories. When the wall tension can't compensate the energy occupied by 
the bubble's volume, it will eventually shrink to zero size. Otherwise the bubble wall will eventually 
expand and quickly approach the speed of light. 

However, these are not the most general solutions. The bubble walls can oscillate and neither 
shrink to zero size nor asymptote to the light cone. This kind of breathing bubble solution, in 
Einstein gravity has been found in the paper [3]. The solution in their case requires the bubble wall 
to have non-standard equation of state. Similar solution is also found in the study of brane-world 
cosmology, where our universe lives on the thin brane embedded in the higher dimension space- 
time. The Israel junction condition here determined the scale factor of the brane-world. In [5], the 
author found that, in order to get the breathing type solution the brane must be embedded in the 
Anti-de Sitter space. 

The Gauss-Bonnet gravity is a particular higher derivative correction to Einstein gravity, which 
leaves the equation of motion remain second order. It is topological in four dimension and will not 
affect the dynamics. In this paper we instead study the bubble dynamics in the five-dimensional(5D) 
Gauss-Bonnet gravity. The energy momentum tensor of the bubble wall is assumed to be the form 
of ordinary scalar field. In many cases the solution is qualitatively the same as in Einstein gravity. 
However in a wide range of parameter space, we find the breathing bubble solution for the true 
vacuum bubble with positive vacuum energy. At first sight, it may seem strange to consider the 
bubble in five-dimensional Gauss-Bonnet gravity. However in the UV complete theory of gravity. It 
is generic to have higher derivative corrections to Einstein gravity. When we consider the evolution 
of the early universe, the effect of higher derivative terms can play an important role. Moreover, 
landscape of vacua in string theory [6J also contains those with effectively space-time dimension 
higher than four. The transition between bubbles with different dimension is also possible [7]. 
There is a possibility that our universe tunnelled from such higher dimensional bubbles. Thus it 
may be too hasty to restrict ourself on 4D Einstein gravity. For example, it needs to be more careful 
when discussing the measure problem, the problem to assign the probability to different vacuum 
bubbles. (See [8] for review). As it usually assumes that the bubble wall approaches the light cone. 

In section two, we will first review the bubble wall dynamics in 5D Einstein gravity. Then we 
turn to the 5D Gauss-Bonnet gravity and find the breathing bubbles. We conclude in section three 
and point out the possible effect of breathing bubble solution on the measure problem. In Appendix 
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A, we discuss the difficulty of the thin wall approximation in general higher derivative gravity other 
than Gauss-Bonnet gravity. However for showing the effects of higher derivative corrections, we 
believe that the Gauss-Bonnet gravity is sufficient. 

2 Bubble dynamics in 5D Gauss-Bonnet gravity 

In this section, we first study the bubble dynamics in 5D Einstein gravity. The solutions will not 
be much different from 4D Einstein gravity if we assume the spherical symmetry. Nevertheless the 
results can then be used to compare with the ones from the 5D Gauss-Bonnet gravity. We then 
turn to bubble dynamics in 5D Gauss-Bonnet gravity. One complication is that the analogy of 
particle dynamics in ID effective potential is not obvious in the Gauss-Bonnet counterpart. Here 
we use the phase space analysis to study the bubble wall trajectories. 

2.1 Vacuum bubble dynamics in 5D Einstein gravity 

We start with the 5D Einstein-Hilbert action 

S=^ j cf>x^[R-2K] (2.1) 

where is the 5D Newton constant, and A is the cosmological constant. Assuming spherical 
symmetry for both the bulk space-time and the bubble wall, we can solve the Einstein equation 
and get the bulk metric solution in the following form 

ds^ = -nn,out{r)dt^ + hnM^y^dr^ + r^d^l (2.2) 

where fin,out{''')^s are the harmonic function for the inside and outside of the bubble. 

Given the inside and outside forms of the metric, the bubble wall dynamics is dictated by the 
Israel junction condition [2], which can be derived from Einstein equation and takes the form as 

[Kab - Khabt = -I^^Sab (2.3) 

where Kab is the extrinsic curvature of the bubble wall, and Sab is the stress tensor of matters on 
the bubble wall. According to spherical symmetry, we can choose the induced metric on the bubble 
wall in the following form 

ds'^ = -dT^ + r^{T)dnl (2.4) 
Then, the stress tensor of the wall takes the form 

Sab = -<yKb (2.5) 
where a can be understood as the tension of the wall. 
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Figure 1: Effective potential for true vacuum bubble in Einstein gravity. 

Choosing the gauge t = t, and plugging the wall metric into the Israel junction condition to 
get the dynamical equation for the bubble radius r{t). Due to the spherical symmetry, the only 
independent equation (the tt component only) takes the form 

- = + f^^^ _ 7^2^ ^ (2.6) 

where "f" '■= This takes the same form as for the 4D Einstein gravity. 

To be more specific, we now consider the dynamics of the true vacuum bubble, that is 

/m = 1 , fout = 1 - - 4 ' (2-7) 

which describe an Minkowski bubble surrounded by the de Sitter space, and ^ is the mass of the 
bubble measured by an asymptotic outside observer. After some rearrangement, (j2.6p can be put 
into form of particle dynamics 

lr2 + K// = (2.8) 

with an effective potential as 

1 , ((4^-A),.^ + Mfrf 
=2- 8 

The form of the effective potential is shown in Figure [H This is quite similar to case for 4D Einstein 
gravity. 

By inspecting the form of the effective potential, it is easy to see there are three kinds of 
solutions of ()2.8p and (j2.9p characterized by the behaviors of r(t). They are bounded, bounce 
and monotonic ones, which correspond to the shrinking, bouncing and ever-expanding bubbles, 
respectively. Again, these solutions are similar to the ones in 4D Einstein gravity. 

One can also change the particle dynamics to the Hamilton formalism by regarding (j2.8p as a 
Hamiltonian constraint, i.e., H = 0. By tuning the initial conditions, one can then numerically 
characterize the above 3 different solutions in the phase space {q = r,p = f) as shown in Figure 
[2jlt seems redundant to numerically solve the bubble wall dynamics in the Hamilton's formalism. 
However, it turns out to be quite essential for the case of Gauss-Bonnet gravity, in which particle 
dynamics is no longer quadratic and is hard to visualize the solution by just inspecting the effective 
potential. 
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Figure 2: True vacuum bubble solutions in phase space for 5D Einstein gravity. A: bounded 
solution. B: monotonic solution. C: bounce solution. 



2.2 Junction condition of 5D Gauss-Bonnet gravity 

We now move to the bubble dynamics in the 5D Gauss-Bonnet gravity by starting with the action 

S=J^j d^x^[R - 2A + aCcB] (2.10) 
where a is the Gauss-Bonnet coupling and the form of Gauss-Bonnet term is 

Cgb = R^- ^R^.uR'''' + R^.upaR''''"' . (2.11) 

Due to the Gauss-Bonnet term, the junction condition for the bubble wall becomes more involved 
as following [3], and it takes the form as 



[Kab - Khab]t + 2a[2,Jab - Jhab + 2PacdbK^'^]t = -K^Sab (2.12) 



where 



and 



Jab = ^{2KKacKl + K.dK^'^Kab - 2KacK'"^Kdb - K^Kab) (2.13) 



Pabcd = Rabcd + 2Rb[c9d]a " 2Ra[cgb]b + R9a[c9d]b ■ (2-14) 



Again assuming spherical symmetry for the bubble wall with the bulk and wall metrics (j2.2p 
and ()2.4p . respectively; then tt component of the junction condition in the t = t gauge can be 
reduced to 



^ ^ - 2a( 2^ ^ '— - 6 1 + i-^)^ I 
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where we have used p.Sp for the stress tensor of the matter on the wah. It is clear that the "particle 
dynamics" for r(t) based on (j2.15p is no longer quadratic, and in general it can be transformed into 
the following form 

X(r)r2 _^ Yij-y^ + Z{r)r^ + V{r) = 

where the exact forms of the functions X{r)^ Y{f)i Z{r) and V{r) depend on fin{r) and fout{f)- 

Therefore, we can no longer use quadratic particle dynamics analogue to characterize the so- 
lutions of r(t), instead it is far easier to use phase space analysis. Let us parameterize the phase 
space as {q = r^p = r), and the Hamiltonian be 

H = X{q)p^ + Y{q)p'^ + Z{q)p^ + V{q) , (2.15) 
then the Hamilton's equations are 



q = 2X{q)p + 4Y{q)p^ +6Z{q)p^ , (2.16) 

dX o dY , dZ d , 

-P = -^P +7r^ ■ 2.17 

oq oq oq oq 

Given fin,outif), we will then numerically solve these equations to characterize the solutions. 

To find the explicit forms of X, Y, Z and V, we assume the spherical symmetry so that bulk 
and wall metrics again take the forms of (j2.2p and (|2.4p . respectively. However, the inside and 
outside harmonic functions solved the field equations of Gauss-Bonnet gravity are the Boulware- 
Deser-Schwarzschild-de Sitter ones [TU], that is. 



/»(r) = 1 + ^(1-^1 + ^) (2.18) 



where A is the cosmological constant inside the bubble; and 



. / N r'^ , / 4aA Sou. 

foutir) = 1 + - + ^ + Tt) (2-19) 

where A is the cosmological constant outside the bubble, and fi is the mass of the bubble. If A is 
zero, this is the true vacuum bubble. On the other hand, the false vacuum bubble has A = 0. 
Given (|2.18p and (j2.19p . the functions X, Y, Z and V appearing in (j2.15p for Hamilton formalism 

are 

t^Cj ly*^ 

maK^i 48aep , QGaA/i ^^ ,2 . , ^^,,0 ^^ ,,00 , , . 2 o„,n;.2.2 






BaA^r^ + SG/cV^ - SaAfcV^ + IGaAAr^ - SaXk^r"^ - SaX^r"^ 

2 



, / o 2a 8 2 .2 , o 2x 8 2 X a 32aA^ 96a^^ 16aA/i\ / 4aA 8a^ 
+ (^_2r A--r aA + 2r A + -r aAA - __ + __j y 1 + __ + __ 

+ + ^ + 2A.^ - 2A.^ + ^aAA.^ - '-aX^r^) , (2.20) 
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, 2a2 ,2 2.^ 2>2 7Q8a^k'^ 192a2A/x ima^Xu STGa^u^ / 



^2 ^4 ^4 

256a^fc^ 

^('") = — ;:2 — ' (2-22) 

-IGa^A^ + 192aA;2 - ^ + 32a2AA - IGa^A^ - 16Au - ^ - 12A;V + 12Au - ^ 
2a 3 

32a//2 288a/i2 192a2A;U ^ ima^X^ A^i? IGaA/x^ _^ 2560^^2 96aA/i 48aA;V 

^6 fy*^ ^4 1^4 ^2 ^2 ^2 

_^96aA/X _ g^^2^2 _^ 3g/.2^2 _ 8qA/.2^2 _^ leo^^jS^^S _ 8q,/.2 -y^2 _ 8q,_)^2^2 _ ^ 

_fi!L _ f a2^4 _ ±„A3r4 + - 2AfcV4 - fcV^ - — + 2AAr4 - 2A:2Ar4 - -W^ - -aX\* 
a 3 27 a 4a 3 27 

,2 8 , , 2 Ar^ 2,24 ^^'^^ 2 , , 2 4 Ar^ 2 , , / 4^^A 8a/7 

+2Ar2 + -aAAr^ AV^ AA;V^ + + - AAr^ ]\ 1 + + — f- 

3 2a 3 2a 3 2a 3 / V 3 

i/'^'" , A\ , ^^/^ , ^6aA^ , OA 2 o\ 2 8 2 822 

+ I-4A/XH 5- H 5 h 2Ar — 2Xr H — aAAr aA r 

V a J.Z 3 3 

Ar^ A;V Ar^ 2,,. 2.^4 2,2 4\ 4^ 



2a 2a 2a 



/Ai 4^ Ar'' Xr^ 2,, / 4c^ I 4aA Sou 



^a 3 Sa^ 6a 6a 
where the re-scaled tension parameter 

k := K^a. (2.24) 

These forms are so comphcated so that we can only rely on the numerical analysis on the phase 
space to characterize various kind of solutions by tuning the parameters and initial conditions. 

2.2.1 Breathing true vacuum bubble 

We first consider the dynamics of true vacuum bubble, which could be relevant to the global picture 
of string landscape. To implement numerical analysis of the Hamilton equations (j2.16p and (j2.17p . 
we set A = l, A = 0, /u = l and a = 0.01 in (|2.2U|) - (|2.23|) . For generic values of k, we will as usual 
have bounded, bounce and monotonic solutions, but for some regime of k we have the breathing 
bubble solution, i.e., the size of the bubble is pulsating/breathing so that the phase space trajectory 
is periodic. The general numerical solutions on phase space are illustrated in Figure [3l 

To examine more closely on the breathing solution, we tune the re-scaled tension parameter k 
for a given some initial condition, and find that there is a critical value of A;, lower than that there 
will exist the breathing solution. The dependence on the initial condition of bubble's phase space 
trajectory is shown in Figure [H 
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Figure 3: True vacuum bubble solutions on phase space for 5D Gauss-Bonnet gravity. A: bounded 
solution. B: monotonic solution. C: bounce solution. D: breathing solution. 




Figure 4: Tuning the re-scaled tension parameter k for breathing bubble. We fix A = 1, n = 1, 
a = 0.01, and the initial condition p = 5 and q = 3.5 at which all the lines intersect in the phase 
space diagram. The solution corresponding to the almost horizontal line has k = 1, and then the 
value of k decreases for the lines arranged counter-clockwise until k = 0.0253 below that there 
exists the breathing solution, i.e., the elliptical shape ones. 
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Figure 5: Phase space diagram for tuning the initial conditions for the true (Left) and the false 
(Right) vacuum bubbles. On the left, we set a = 0.01, k = 0.027, /i = 1, A = 1 and A = 0. On the 
right, we set a = 0.01, k = 0.001, /x = 1, A = and A = 1. We see there is a breathing solution for 
the true vacuum bubble, but not for the false one. 

2.2.2 No breathing false vacuum bubble 

We can also consider the dynamics of the false vacuum by choosing A = and A = 1 for the 
harmonic functions (|2.18|) and (j2.19|) . Again, solving the Hamilton equations (j2.16|) - (|2.17|) with the 
set of functions in (12.20p - (|2.23p in terms of the new harmonic functions, we find that there is no 
breathing solution. For illustration, in Figure [5] we juxtapose the phase space diagrams by tuning 
the initial conditions for the true and false vacuum bubbles. 

We can also consider the cases with both A and A non-zero, then it is more close to true vacuum 
bubble if A > A, and more close to false vacuum bubble if A < A. From our numerical studies with 
quite extensive ranges for both parameters and initial conditions, we find that there exist breathing 
bubble solutions for the former case, but not for the latter. We may then tend to state that there 
are breathing vacuum bubble solutions in 5D Gauss-Bonnet gravity if the cosmological constant 
inside the bubble is smaller than the one outside, but not vice versa. 

3 Conclusion and the implications 

The higher derivative terms in Einstein-Hilbert action can be ignored when the length scale con- 
sidered is much smaller than the curvature radius. Thus there is no reason to dwell on the Einstein 
gravity when we consider the global picture of the universe. In this case it will be important to 
exam the effect of the higher derivative corrections. One new feature we found in this paper is 
the breathing bubble solution in 5D Gauss-Bonnet gravity. In particular, we don't need the exotic 
matter for the bubble wall and this kind of solutions exist for true vacuum bubble with zero or 
positive vacuum energy. One possible implication is on the measure problem of eternal inflation. In 
order to do the sensible statistics on bubble distribution, we need to cut-off the infinite space-time 
and infinite number of bubbles properly. One way of doing this is to count the bubbles smaller 
than a fixed scale factor [llj. However usually people assume the bubble nucleating at vacuum 
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with Hubble constant H will have co-moving radius H asymptotically. This is not the case 
if there are breathing bubbles whose co-moving sizes are vary with time and depend on critical 
bubble size. Including these bubbles will change the evolution equation for the fraction of volume 
occupied by a particular bubble. One may also ask what happens if a growing bubble nucleates 
inside a breathing bubble. The bubble walls can collide and generate the gravitational wave. This 
can have the observational consequences if our universe goes through such kind of bubbles in the 
past evolution. We leave these problems for the future studies. 
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A Difficulty for thin-wall vacuum bubbles in higher derivative 
gravity 

One difficulty of studying bubble dynamics in arbitrary higher derivative gravity is due to the 
appearance of higher order singularity when assuming the thin wall approximation. In [12], they 
try to avoid this difficulty by assuming the discontinuity of extrinsic curvature only appears in higher 
order derivative. However we don't have solutions satisfing this criterion yet. In this appendix, we 
will derive the junction conditions modified by higher order corrections and argue that in general 
they are not the consistent sets of equations. 

A Lagrangian including higher order corrections derived from high energy theory like string 
theory may have the form 

L = R-2A + aR^ + fiR^^uR^" + jRa/s-ysR"^"'^ 

A special case for second order gravity isa = 7 = l,/3 = —4 which is the Einstein-Gauss-Bonnet 
gravity. In four dimensions these terms are topological invariant which does not contribute to the 
equation of motion. For simplicity, we consider L^^^ = R — 2A + aR^ first and show that thin wall 
approximation is not reliable. The field equation from variation of ^/gL is |12] 

(7^^ = (1 + 2aR)Gf,^ + ^aR^gf,^ - 2aDf,D^R + 2ag^yUR^ 

Uij = Ahij + Gij + 2aRGij + a^R^hij + 2a[-DiDjR + KijdyR + hij{dyyR + Di? - KdyR)] 

In order to extract the discontinuity from equation of motion, we assume the metric is continuous 
at the wall, but it has a kink. Its first derivative has a step function discontinuity and its second 
derivative has a delta function term [TH] 

h^J{y) = h-^ei-y) + h±e{y) 
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In Einstein's gravity, substitute (24) into Gij = K^Tij and match the delta function term, so that 
we can get the Israel junction condition. In order to get the junction condition for the correction 
gravity, we need to work out following quantity 

RG^j = 2dyKdy{Kij - Khij) + 2dyK{2K\Kij + dyKij - KKij + ^KjiTrK^ + K^) + G,,) 
+ dy{Kij - Khij){-TrK^ -K^ + R) (A.l) 

R^ = [2dyK - TrK'^ - K'^ + R]^ (A.2) 

dyRKij = [2dyyK - dyTrK^ - dyK"^ + dy{R)]Kij (A.3) 

dyyR = 2dyyyK - OyyTrK^ - OyyK^ + dyy{R) (A.4) 

KdyR = K{2dyyK - dyTrK"^ - dyK^ + dy{R) (A.5) 

We consider that induced metric hij is continuous across the domain wall, but discontinuous across 
the wall. However terms like R^ would give 5{y)'^ and dyR would give dy6{y) and dyyR would give 
dyy6{y). Those more singular terms are not well defined. One may require higher continuity of 
the induced metric hij, but this would make the equation of the domain wall trivial. Wc cannot 
assume this condition in our case. We propose to include more singular source terms on the wall 
to see if the equation of motion can be solved consistently. Introducing source terms 

= S^''{x')6{y) + U^''{x')6{yf + V''''{x')dy6{y) + W^'''{x')dyy6{y) 
Due to the conservation lawVi^T'*'^ = 

V^T'" = {VjS'^ + 2KjS^y + TrKS^)5{y) + S^dy5{y) 
+ 2K'jWy + TrKU^y)5{yf + U'y2d{y)dyS{y) 
+{VjV'^ + 2K'jV^y + TrKV'y)dy5{y) + V'ydyy6{y) 
+{VjW'^ + 2K'jW^y + TrKW'y)dyyS{y) + W'ydyyy6{y) 

= 

We should have the following constraints. 

VjS'^ + 2Kp^y + TrKS'y = 
VjU'^ + 2Kp^y + TrKU'y = 
S'y + VjV'^ + 2KiV^y + TrKV'y = 

u^y = Q 

VjW'^ + 2Kiw^y + TrKW'y + v'y = Q 

w'y = Q 
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Matching the left hand side and right hand side of equation for the first and second derivatives of 
the delta function and the square of the delta function respectively. It gives four equations, where 
one being proportional to 6{y) is rather cumbersome, so we neglect it, but it will not affect our 
result 

dyySiy) -2ah'^+{dht - dhT.)hij = Wi,{x') 



\[h''^{dyKl - dyKi)h+'^''idyh+^ " dyh-Jidyh± + dyh;^)] + ^[h^'' {dyh+ - dyhT.)fh,, = U., 



dyS{y) 2^h+''dyKl^dyhl, + h+^'dyhT^dy + h+^^h^^' dyh+^dyhl, + h^^^ dyhT^dy + h,,)h,, 

- 2{dyh + ^^dyh+ + dyh^^^dyhr^ + h-^'^ d^yh^ + h-^'^ dyhT. + 9, /i^^^' (9, /i+ - dyh-.) 
+ h-^'^idyh+ - dyh-^) + dyh~^^^dyh± + h-^'^d^hf. - dyh'^^ dyK^ - h-^^ dyh-^)h,J 
+ ih^^'h+^\dyh±dyhtl - dyh.^dyh-) + h^^^ h + '^dyhf^dyht, - dyh-^dyh~,))h.,] = V^J 

We can see that these three equations are very different and cannot be solved simultaneously. It 
shows the breakdown of the thin wall approximation. Instead, we should consider finite thickness 
of the domain wall to avoid severe singular behavior. 
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